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velocity fields of a black hole and a rigid sphere in an external dipolar field (constant acceleration 
field). In the latter case the dipole field is produced by an axially symmetric halo or shell of 
matter. For each case the fluid density is studied using contour lines. We found that the presence 
of acceleration is detected by these contour lines. As far as we know this is the first time that 
the integral curves of the velocity field for accelerate objects and related spacetimes are studied in 
general relativity. 
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I. INTRODUCTION 



The study of potential flows in general relativity is relevant in the understanding of several phenomena of interest 
in relativistic astrophysics like: fluid flows at relativistic speeds in the presence of a neutron star [1,2], star clusters 
moving in gaseous media [1], flows near a cosmic string [1], accretion in binary star systems and supermassive black 
holes [3], and others [4-6]. Also, numerical solutions of the equations of general relativistic hydrodynamics can 
■ simulate, and model, gravitational collapse and the evolution of neutron stars [7]. 
t-H | Most of the articles in this area deal with the important case of fluid motion evolving in the spacetime associated 
with compact stars and black holes. The implementation of new background metrics brings some new challenges. 
"q I First, metrics other than Schwarzschild and Kerr are not so well studied, sometimes a complete understanding of the 
physical meaning of the metric is missing. Also, the solutions of the fluid equations in a non trivial metric may be 
quite involved. In particular, the search for significant boundary values (or initial conditions) presents a non trivial 
j^Tj, problem. A simple and paradigmatic case of a flow is the stationary-zero-vorticity flow of a fluid with adiabatic stiff 
equation of state. In this case, for relativistic flows, the fluid equations admit analytical solutions for some particular 
metrics [1,3]. These solutions are used as test-beds for testing almost all the numerical hydrodynamic codes in the 
subject. Other potential flows in a non stationary background and different equation of state have been studied, 
see for instance [2]. Also, the solutions for potential flows permit to test new optimized codes in resolving nonlinear 
hyperbolic systems of conservation laws [6] , see for a representative sample of numerical schemes [7] . 

In this work we extend the investigations about potential flows by studying the streamlines of a steady-state ideal 
fluid in the presence of an accelerated black hole and rigid sphere, and of a black hole and rigid sphere with a dipole 
field produced by an axially symmetric halo of matter. This shell like structures are useful in modeling many situations 
of interest in astrophysics, as for example, the Supernova 1987A [8-10], for other applications see [11]. We assume 
that the fluid is a test fluid, i.e., the metric does not evolve and it is given apriori. The state equation and the idea 
of a rigid star (rigid sphere) we use are idealized. However, they bring important results about the instability and 
behaviour of the fluid, and the difficulties involved in this kind of scenarios. 

This work is divided as follows. In Sec. II we present the basic equations that describe potential flows. In Sec. 
Ill we summarize some aspects of the Weyl C-metric that represents the spacetime associated with an uniformly 
accelerated black hole. In particular, we present the metric in different systems of coordinates to facilitate the 
physical interpretation of the results. In the subsections III. A and III.B we study potential flows for an accelerated 
black hole and an accelerated rigid sphere, respectively. In the first case we use a perturbative approach and in the 
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second we solve numerically the potential equation. For both cases, we also study the behaviour of the fluid density. 
In Sec. IV we introduce a metric that represents a black hole in a dipolar field, this field is produced by an external 
halo of matter. In the subsections IV. A and IV. B, respectively, we study the streamlines of a fluid in the presence of 
a rigid sphere and black hole, both with halo. Also, in both cases we study the fluid density. Finally in Sec. V we 
summarized our results. 



II. BASIC EQUATIONS 



The starting point of this work is an ideal fluid whose energy-moment tensor is given by = {p + pjU^U,, + pg^ v , 
where p is the pressure, p the total energy density and the four-velocity. The conservation equations, T.£" = 0, for 
this kind of fluid reduce to 



and 



{p+p)m + p, li u'* = o, 



{p + pWU^+p^+p^WU^ 



0, 



(1) 

(2) 



which are respectively the conservation an Euler equation. (Our conventions are: G = c = 1. Metric with signature 
+2. Partial and covariant derivatives with respect to the coordinate x M denoted by , p and ;p, respectively.) For 
isentropic flows we have, (a/n)-^ — 0, where a is the entropy per unit volume and n the baryon number density. In 
this case the equations of motion (2) take the form [12], 



0. 



where is the relativistic vorticity tensor defined as 



l^uv — 



[("?K 













(3) 



(4) 



The potential flow solution of this equation {w^ v = 0) is 

' p + p 



where $ is a scalar field. From (5) and the equation of continuity for the baryon number density n, (nU^).^ 
obtain the differential equation for the scalar field $, 



p + p 



0. 



(5) 
0, we 

(6) 



The normalization condition, U^U^ = —1, provides a relation between the pressure, the total energy density, and 
scalar field: (p + p)/n = < E > ,^ ( E >,M - Also, p and p are assumed to be related by a barotropic equation of state, 
P = p(p)- I n general, Eq. (6) is a nonlinear equation for the scalar field However it becomes linear if we assume 
that p = p oc n 2 , i.e., a stiff equation of state. Thus, 



nU^ = $ M , 



n 2 = CT $ ,<T , 



and $ is a solution of 



□* = (>/=0<r $,„),„ = o, 



(7) 

(8) 



which is the usual wave equation for a scalar field. Note that for the stiff equation of state the sound velocity in the 
fluid is equal to the velocity of light, therefore the flow is always subsonic and the presence of Shockwaves is excluded. 



2 



III. ACCELERATED BLACK HOLES AND RIGID SPHERES 



The Weyl C-metric is the member of the static axially symmetric Weyl family of solution of the vacuum Einstein 
equations [13], 



ds 2 



-F(q)dt 2 + 



dp 2 



dq 2 



G(p) F(q) 



(9) 



(10) 



A 2 (q+p) 2 

where the functions G{p) and F(q) are the cubic polynomials, 

F(q) = -1 + q 2 - 2mAq z , 
G(p) = 1 -p 2 - 2m Ap 3 . 

The coordinates t, p, q and <p are dimensionless, A is a constant with dimension of L~ x . The range of t is (— oo, +oo), 
while the range of ip is [0, 2ir]. The ranges of the coordinates p and q depend on the roots of F and G. The constraint 
m 2 A 2 < 1/27 is imposed to ensure the existence of three real roots in (10). This condition is not a physical one, it is 
only due to the choice of coordinates [14]. From the fact that the functions have three real roots, the C-metric allows 
a description of different spacetimes [14] depending on the interval considered for (q,p). In this work the functions 
F(q) and G(p) are bounded within the zeros (qn,qs) and (po,Pn) given in [15]. In this case , m and A are interpreted 
as the mass and acceleration of a black hole [16]. A Newtonian image of the matter content of this metric is given 
by a rod and a semi-infinite line mass, both of density 1/2 and placed along the z-axes. The metric is interpreted as 
the representation of particle (the rod) in a accelerating frame (the semi-infinite line), in other words, associated to 
the spacetimes of Schwarzschild and Rindler. The rod and the semi-infinite line has a strut holding them apart that 
represents a conical singularity. We note that the metric (9) diverges when the acceleration is zero (A —> 0). Hence, 
in this coordinates, we do not have the correct limit that corresponds to the Schwarzschild metric. To obtain this last 
metric when A ^ we make the coordinate transformation [15], 

1 

r = — -, t-fAt. 



Now the metric (9) reads, 



ds 1 



A{q+p) 



-Hdt 2 + \zdr 2 + ^—drdp + r 2 
H H 



dp 2 + r 2 Gdip 2 



(11) 



(12) 



where 



H = A 2 r 2 F 

= -A 2 r 2 G(p-A- 



= 1 - 



2m 



+ 6Amp + ArG^p - A 2 r 2 G(p) 



(13) 



When A — > the metric (12) reduces to Schwarzschild form provided that the spherical angular coordinate 9 be 
related to p by G{p) = 1 — p 2 = sin 2 9. Furthermore, if the mass is zero in (12) ( m — > 0) the space becomes Euclidean 
with a special form of a flat space in an uniformly accelerated frame [15]. For these reasons the line element (12) 
represents a uniformly accelerating Schwarzschild-type particle. The coordinates (i, r, p, if) is a coordinate system 
rigidly fixed on the accelerating particle. With the transformation (11) we gain the correct Schwarzschild limit, but 
now the metric is not diagonal. 

In the case of A ^ 0, the angular coordinates p and 9 are related by G(p) — sin 2 9. The mapping between them is 
[15], 



, = (-_ L 

\ 6Am 



(Q{6) 4tt 
2cos ( 3 + T 



for < 9 < 



6Am 



2 cos 



0(9) 



2tt 

Y 



+ i 



< IT, 



where 



cos 0(9) = 1 - 54A 2 m 2 cos 2 9. 
Due to the acceleration the black hole horizon deforms. The exact form of this horizon is [15], 



(14) 



(15) 



3m 



rsch 



cos 



+ 



4 7l 



)+cos(f + f) 



forO < 9 < -, 



3m 



cos (®M + 



¥)+cos(| + f)" " 



for - < 9 < it. (16) 



The form of the Schwarzschild horizon will be used in a perturbation scheme to place an upper limit on the acceleration. 



3 



A. Potential flows and accelerated black holes 



In order to study streamlines of the flow we need to know the field $ that determines the fluid velocity. The explicit 
form of the equation for <J> can be obtained from (8) and (12). We find, 

□$ = ~$ tU + ({H + A 2 r 2 G}rH. r ), r + (G$. p ), p + - (r 2 GA<S>. r ). p - (r 2 GA$ p ), r = 0. (17) 

rl (_x 

We shall assume that the fluid is stationary and also independent of the variable ip. In this case the time dependence 
of $ will be set by adding the term — at to the final solution, where a is a constant, this term clearly satisfies (17). 
The arbitrary constant a is related to the th component of the four-velocity. With these assumptions, and setting 
m = 1, (17) takes the form 



with 



□$ = a Sch $ + Aa 1 ^ + A 2 a 2 ^ = o, (18) 
□sch* - ([r 2 - 2r]$, r ), r + ([1 - p 2 ]$ P ),p, 

□ i$ = ([6pr 2 - 2pr 3 ]$ r ), r - (2p 3 $ p ), p - (r 2 [l - p 2 }$, r ), P - (r 2 [l - p 2 ]$, P ), r , (19) 

□ 2 $ = -(6rV*,r),r + (2rV*,r)j, + (2r 2 p 3 *j.),r- 

Finding an analytical solution to Eq. (18) does not appear possible. We shall look for a meaningful approximate 
solution. For small enough A, say A < 0.01, the Schwarzschild surface remains almost unaltered [cf. (16)]. Then we 
can approximate (18) as 

□ $ = n Sch <$> + ADi$ = 0. (20) 

By using the fact that A ss we may solve this equation by a perturbative method, but first we recall some earlier 
results. The solution for a black hole when A = with the condition that at infinity the fluid velocity is constant and 
parallel to the z-axes of the inertial frame is known [3] , 

$ PST = -at - 2a In ^1 - ^ + b(r - l)cos0, (21) 

where a and b are constants. We shall consider for the unperturbed potential (A = 0) that the fluid is at rest at 
infinity. Following the procedure of [3] we find that the unperturbed potential is, 

$o = -at - 2a In ^1 - . (22) 

This solution can also be found in a different way. Note that $. r diverges near the black hole horizon. This is a 
necessary condition to avoid the divergence of the particle density on the horizon [3]. Therefore, near the black hole 
horizon, we can neglect <& iP when compared with It is illustrative to consider (18) with A = 

□scfc* = ([r 2 - 2r]$ r ), r + ([1 - p 2 ]$ P ), P = 0. (23) 

When near the black hole we neglect the terms containing derivatives with respect to p the equation reduces to an 
equation for r only. By considering that the temporal dependence is given by —at this equation has for solution (22). 
Hence, the above imposed condition is valid and moreover, in the case of fluid at rest at infinity, gives the exact 
solution. 

Now we shall consider the perturbation, 

$ = $o + ^4T, (24) 
where T is a function of r and p. From (20) we find, 

a Sch T = -n&o. (25) 
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Thus, 



([r 2 -2r]T ir ) ir . + ([l-p 2 ]T, p ), p 



16ap(r — 3) 
(r-2) 2 ' 



(26) 



Again, to solve this equation we use the fact that near the black hole we can neglect the term containing derivatives 
of p. We find, 



t (■ ■ ir,:.: (\ H ) _ i up 



2 , (■> 2 

o + In 1 - - 

r - 2 \ r 



(27) 



+2 In (0 ln(r - 2) - [ln(r - 2)] 2 + 2PolyLog (2, 1 



with 



PolyLog(n, z) = ^ 



oo t. 
Z K 



fc=l 



(28) 



To obtain the range of validity of the perturbation AT, we compare it with the unperturbed function <J>o- We 
consider that the perturbation is valid only when the absolute value of AT is at least a 15% of the value of $o- We 
can take values of the constants G\ and C 2 to enlarge the range of the radial coordinate and still be in the required 
precision. For example, for A k> 0.01 and the G\ and Ci in the interval (0, 2) the range of the radial coordinate is 
r sa (2, 10). In Fig. 1. we show the flux lines for the case A = (dotted lines) and A ss 0.01 (full lines). When 
A w 0.01, the streamlines that were straight lines are now curved due to the black hole acceleration. These lines are 
similar to the corresponding lines for a fluid flow in the presence of a black hole, with the condition that the fluid 
velocity be constant at infinity [1,3]. 



Another significant quantity is the particle density of the fluid, n 2 = — At first order in A we obtain 



16a 2 



n 



1 + 6Ap 



+ 



8aAC 2 



2Apr (r - 2)r 3 (r - 2)r 



+ 32aVl 



r(r - 3) + 4{1 + (r - 2) ln(r - 2)} 



(r-2) 



(29) 



The particle density (29) is positive in the interval wherein the perturbation is valid. When A — > we obtain the 
correct limit for a fluid at rest at infinity [3]. The fluid density contours are shown in Fig. 2. We see that near the 
black hole these lines are closed and more dense in the forward direction. These contour lines are also similar to the 
corresponding density lines of a fluid in the presence of a black hole, with the condition that the fluid velocity be 
constant at infinity [3]. Far from the black hole these lines are open, while in [3] they are closed, this clearly shows 
the difference between a standing black hole in a moving fluid and an accelerated black hole moving in a fluid which 
is at rest at infinity. The qualitative aspects of the density contours do not vary when lowering the acceleration or 
varying the constant C2. When A — > we density contours are circles centered in the black hole. 



B. Potential flows and accelerated rigid spheres 



In the case of rigid spheres moving in a fluid, the particle density is no longer divergent on the sphere surface nor 
<f> r is divergent. Then, even near the surface of the sphere, we cannot neglect terms containing derivatives of p in the 
equation for <f>. To find a numerical solution we need to solve the partial differential Eq. (8) with mixed boundary 
conditions, Neumann and Dirichlet. In this case it is more convenient to work with the metric (9) which is diagonal 
and later perform the coordinate transformation (11). In this metric Eq. (8) reads, 



F{q + pf 



G 

(q + p) 2 ' P . 



F 

(q + p) 2 '". 



G{q+pf 



0. 



(30) 



From the condition that the fluid is stationary and that $ is independent of ip we get, 

F 



G 

(q + p) 2 ' P . 



= 0. 



(31) 
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The boundary conditions are: i) Zero normal component of the fluid velocity on the sphere surface, and ii) Faraway 
from the sphere we set the fluid velocity to be constant and parallel to the acceleration. The condition i) do not 
describe the typical flow of gas around a star because the surface of the star is not hard but gaseous, except in special 
astrophysical situations; but describes an idealized strong-field star. 

The code used to solve (31) is a finite difference multigrid method [17] with second order precision. The numerical 
multigrid is evenly spaced in r, but not in 9. We find that the convergence is better when we increase the number 
of points in 9 rather than in r. To solve (31) we employ four domains (grids) limited by spherical shells of radius 
greater than 2.5 (Schwarzschild radius equal 2). For the first grid we have r <G [2.5, 3.9] and we use 79300 points; the 
respective numbers for the other three grids are [3.9, 5.3], 39700, [5.3, 11.1], 19900, and [11.1, 25.5], 5000. An iterative 
method in which at every step each point of the grid is calculated from the values of the four nearest neighbors is 
employed. The error is estimated with the expression \&new — ^oid\> m which $^ is the old value of $ at the point 
{xi) and ^nel is the new calculated value at the same point. The program will stop when the sum of all errors of 
the grid points reach some pre-establish value, say ^2 X . l^n^'i — ^Im^I < Error. The code was tested with the exact 
solution for a rigid sphere in a fluid flow [1]. The numerical results on the surface of the sphere have less than 1% of 
error. Since we are interested in qualitative aspects of the fluid behavior rather than in precise figures this accuracy 
is sufficient for our purposes. The precision in the rest of the grid is better. 

In Fig. 3. the streamlines of a fluid disturbed by a moving rigid sphere with constant acceleration A = 0.01 are 
depicted. We see that the streamlines concentrate in the frontal part and separate when passing the rigid sphere. This 
may be due to the deformation that suffer the killing horizons in an accelerated frame [15]. The qualitative aspects 
of the density contours for this case are similar to the black hole case. 



IV. POTENTIAL FLOWS IN BLACK HOLES AND RIGID SPHERES WITH DIPOLAR HALO 



To incorporate the dipolar field into the Schwarzschild metric we consider the static axially symmetric Weyl metric, 

ds 2 = -e 2 »dt 2 + e 2 ^ (dz 2 + dp 2 ) + e- 2 »p 2 d V 2 , (32) 
where v and 7 are functions of p and z only and satisfy the conditions [18], 

1 



v.pp + -v p + v zz = 0, 
p 



dj = p [(^p) 2 - {v,z) 2 } dp + 2pis iP v tZ dz. 



(33) 
(34) 



The first equation is the usual Laplace's equation in cylindrical coordinates, and the second (once v is known) gives 
7 as a quadrature. 

For a multipolar expansion of v 1 the spherical coordinates (r, 9, ip) or the prolate spherical coordinates (u, v, (p) are 
more adequate than the cylindrical ones. The relation between these coordinates are 

U = 2ro Vp 2 + ( z + m ) 2 + V P 2 + ( z - m ) 2 
r 

= 1, u > 1, 

TO 
1 



2to 
cos 9 



V 'p 2 + (z + to) 2 - V 'p 2 + {z- m) 2 
-1< v < 1, 



(35) 



ip = <p. 



In terms of p and z, 



p = m\J (u 2 — 1)(1 — v 2 ) 
= \J r(r — 2m) sin 9, r > 2m, 

z = muv 
= (r — to) cos 9, 

(p = <p. 



(36) 



As before, we will set to = 1. 



G 



Using the transformations (36), Eqs. (33) and (34) can be written in terms of u and v as 

[(« 2 -1K«] iU +[(1-«>,v] i „ = 0, (37) 

7,u = -5 5- [ M (w 2 - l)( I/ ,«) 2 - M (l - w 2 )(^.d) 2 - 2w(u 2 - l)i>, u i>, v ] , 

u v 

u 2 - 1 

7,« = -5 o [«(« 2 - 1)0.«) 2 - «(1 - ^X".*) 2 + 2u(l - w 2 )i/, u i/,v] ■ 

u v 

The authors in [11] solved this equations using a external multipolar expansion up to octopoles using a Legendre 
expansion with their corresponding terms increasing with the distance in the intermediate vacuum between the core 
and the shell. In this work we are only interested in the dipolar case. Hence, the functions v and 7 reduce to 

2u = 2u + «ln ( — - ) + 2Vuv, (38) 
\u+lj 

2 7 = 2 7o + k 2 In ( JL-ZL ^J + AnVv - V 2 [u 2 {l - v 2 ) + v 2 ] , (39) 

where T> is the value of the dipolar field produced by an axially symmetric halo (or shell) of matter, and n is a 
constant. The Schwarzschild solution is recovered with k = 1 and T> = 0. The constants 70 and vq can be used to 
rule out conical singularities and to ensure analyticity of the metric at the horizon [11]; here we have made them zero. 
From the metric (32) in coordinates u and v, and (38) and (39), we find 



ds z = - 

+ 



" l\ e 2T> uv dt 2 + ( u+ l \2 e 2Vv(2-u)-V 2 [u 2 (l-v 2 )+v 2 ] 

U+lJ 



U + 

dv 2 



du 2 



l-v 2 



u 2 - 1 

+ (u+l) 2 (l-v 2 )e- 2Vuv dcp 2 . (40) 



This metric represents a monopolar core (k = 1) in the presence of a external dipolar field (V) that is associated to 
a distant shell or halo of matter. 



A. Potential flow of a rigid sphere in a external dipolar field 

For a stationary fluid that does not depend on ip the scalar equation (8) reduces to 

[(« 2 -l)*, u ] itt +[(l-t; 2 )$,„] ii; = 0. (41) 

This is exactly the same equation that we solved in Sec. III. A, but in a different system of coordinates. Note that all 
the Weyl solutions have the same differential equation (41) for $. The presence of the dipolar field will be taken into 
account in the boundary conditions. On the surface of the rigid sphere the boundary condition is the usual one for 
the Euler's equation, i.e., that the normal component of the fluid velocity vanish on the surface. For large values of 
u the boundary condition deserves more attention. Since, the time dependence of <& is —at, we have for the density, 

/„ 1 1\ 2Vv(u-2)+V 2 [u 2 (l-v 2 )+v 2 ] 

n 2 = a 2 (!l±lj e~«^ - ^—y 2 [(«,„) V - 1) + (*,„) a (l - « 2 )] • (42) 

Due to the presence of the dipolar field T>, at some value of u greater than a certain u = u c , the first term in (42) 
will always be smaller than the others two. Then we will have n 2 < 0. This is not allowed; we say that the fluid is 
no longer stationary, i.e., the assumed time dependence is not right for u > u c . If T> is large enough, say Pw 1, the 
fluid will never be stationary. The first term also depends on a 2 . By varying the value of a we can enlarge or decrease 
the domain where the fluid is stationary, also for V — > 00 we can set a large enough to keep the fluid stationary. For 
example, for a = 1.25, we must have T> w 0.001 to obtain a relative large domain wherein the fluid remains stationary. 
Test particles moving in the metric (40) can have chaotic behaviour [11]. We believe that the lost of the stationary 
character of the fluid for u > u c may be another manifestation of the same type of instability. This point is under 
active consideration by the authors. 

The outer boundary condition is still missing. With all these fluid instabilities it is not clear which is the right 
boundary condition. The previous condition, that the fluid velocity to be constant far from the sphere, is not a valid 
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condition in this case because the fluid is accelerating and is in this region where the instabilities appear. To find some 
characteristics values of $ to be used as an outer boundary condition we will integrate $ jtl along a line of constant v 
from the surface of the sphere (black hole) to the stationary limit of the fluid. From Eq. (7) we get 

(43) 

The proportionality constant is of no importance, because in (7) is invariant under a rescaling of $. To find the 
t component of the four- velocity, [/*, we will assume - without further justification - that the a fluid particle follows 
the usual geodesic equation of motion for a test particle, i.e., the Euler-Lagrange equation, 

where x^ represents total derivative of x^ with respect to the parameter s and L(x x ,x x ) — ^g^x^x". Since the 
coordinate t is cyclic, we find, 



t = U = —k y^—^J e , (45) 

where the constant fc is the energy function of the fluid particle. The component U u along the lines of constant v can 
be computed from the metric (40), 



U u = u = J( ) e^uv + fc 2 e-mv+^W-v^+v*] ( 46 ) 



,u+ 1 

Therefore, $ along the lines of constant v can be written as 




where a multiplicative constant was set equal to one. We numerically integrate (47) to obtain the outer boundary 
condition at a certain value u near the limit of stability. We use the code described in Sec. III.B to solve (41). For 
different values of the constant k the qualitative features of the streamlines and density contours do not change. In 
Fig. 4 the streamlines of a fluid in the presence of a rigid sphere of radius 2.5, Schwarzschild radius equal 2, and a 
external dipolar moment field with T> = 0.001 are shown. We also set fc = 1 in this case. The values of $ used for 
boundary condition and all the ones computed satisfy the condition n 2 > 0. We see little difference in the flows with 
T> = 0.001 and T> = 0. In Fig. 5 the density contours are plotted for the same values of the parameters used in Fig 4, 
sphere of radius 2.5 and V = 0.001. We see the same qualitative features discussed in Sec. II. A in the context of the 
C-metric. In this case the fluid is accelerating something that we should expect since an external field (in this case 
the external dipolar field) exert a force on the fluid. 



B. Potential flows for black holes with dipolar halo 



Like in the precedent case we solve Eq. (41), but to avoid the singularity at the black hole horizon, we use the 
tortoise radial coordinate [19], 

r* =r + 21n(^ - l) (48) 

= (« + l) + 21n(^), (49) 

Now, Eq. (41) reads, 

[{u + 1)2 * r * ] + [(1 - w2)$ - ] .« = °' (50) 

and the particle density takes the form, 
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2 2 

n = a 



1 



-2Vuv 2Vv(u~2)+V 2 [u 2 (l-v 2 )+v 2 



" + 1 -2 1 — 



It- 1 



($, r .) + 



(51) 



From a computational view point the difference between a rigid sphere and a black hole laid in the inner boundary 
condition. For a rigid sphere the condition is zero normal velocity on its surface, and for a black hole is finite particle 
density on the horizon. This last condition requires that <& ir * be limited in such a way that cancels the singular term 
in $ „. In the general time dependent case this condition reads [2,6], 



1 



Hi * 

r , 



$, t + oi(t)(r - 2) + a 2 (t)(r - 2) 2 + 



where a\, a 2 , ... are functions of time coordinate. Hence, the following equation is valid near r = 2 

0. 



_d_ 

dr* 



r - 2 



(52) 



(53) 



With the assumption that <& it = —a this condition is also valid in our case and it will be taken as the inner boundary 
condition. In (52) we cannot set <&. r . = <& it because we can have n 2 < 0. The outer condition for this problem is 
found integrating Eq. (47) as before. The numerical code employed is the same of Sec. IV. A with the implementation 
of the inner boundary condition (53) and the change of r or u to the tortoise coordinate r* . The program was tested 
with the exact solution for a black hole in a moving fluid [3] . The computed values of <& on the black hole horizon 
have less than 1% of error when compared with the exact solution. Again, outside the horizon the precision is better. 

In Fig. 6 the streamlines of a perfect fluid in the presence of a black hole in an external dipolar moment field with 
V = 0.001 are plotted. Like the rigid sphere case the streamlines show little difference with respect to the case V = 0. 
In Fig. 7. we present the density contours for the same values of the parameters used in the previous figure. Once 
more we see the same pattern encountered in the C-metric case. This shows that the effect of acceleration appears 
always in the density contours rather than in the fluid streamlines. In the last four figures we see that the streamlines 
are practically the same compared to the respective case with T> = 0, but the density contours are quite different. 
The density contours for the different obstacles present remarkable similar features. We think that this pattern may 
be taken as test for the presence of acceleration in this kind of fluids. 



V. CONCLUSIONS 



We can summarized the results of this work as follows: in the C-metric and dipolar field cases we found that, for a 
stationary fluid, the difference between the accelerate and non accelerate cases lies on the form of the density contours 
and not, as one may think, in the shape of the streamlines. We assumed in both cases that the acceleration is near 
zero. In the dipolar field case the potential flow becomes very unstable and the presence of chaos may appear [11]. 

It is important to note that the method presented hear is valid for a stiff equation of state and the incorporation 
of a new barotropic equation of state for the fluid is not easily applied. In that case we need a different approach 
for solving the difference equation (6) because it is nonlinear and it is necessary to impose initial conditions to the 
enthalpy, this may help or not in the convergence of the method. A negative value of the enthalpy could indicate a 
nonstationary regime for the fluid , see for instance [2] . In the different scenarios of this work the search for suitable 
boundary conditions to be treated, with 4«0, have been a constant problem, we think that this difficulties would 
increase with a new barotropic equation due to the condition for the enthalpy. Future applications of this work could 
be the study of fluids with another equation of state to model a more realistic situation, the study of chaos in the 
presence of instabilities, and the study of fluids in non stationary states. All this applications are currently under 
investigation by the authors. 
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FIG. 1. Analytical results for the streamlines of a fluid in the presence of an accelerating black hole in the direction of the 
positive 2-axis. The solid lines represent the case with A 7^ and the dotted lines A — 0. The black hole has radius r ~ 2 
(Schwarzschild radius equal 2), this approximation is due to the deformation of the Schwarzschild horizon. The axes are defined 
as X — r sin 6 and Z = - - - •■. ... _ , 




FIG. 2. Analytical results of the fluid density contours of an accelerating black hole. The black hole radius and the meaning 
of the axes are defined as in Fig. 1. 
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FIG. 3. Numerical results for the fluid streamlines in the presence of a rigid sphere accelerating in the direction of the 
positive z-axis. The sphere radius is r ~ 2.5 (Schwarzschild radius equal 2), this approximation is due to the deformation of 
the Schwarzschild hori 




FIG. 4. Numerical results for the fluid streamlines when an external dipolar field with T> = 0.001 is present and a rigid sphere 
of radius r — 2.5 (Schwarzschild radius equal 2) is placed as an obstacle. The axes are defined as X = rsin^ and Z = rcos#, 
with r = u + 1 and cos 6 = v. 
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FIG. 5. Numerical results for the fluid density contours when an external dipolar field with T> = 0.001 is present and a rigid 
sphere of radius r = 2,~ 




-i& X is 



FIG. 6. Numerical results for the streamlines when an external dipolar field of value T> = 0.001 is present and a black hole 
is placed as an obstacle. The black hole has radius r = 2 (Schwarzschild radius equal 2). The meaning of axes are the same of 
Fig. 4. 
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